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ABSTRACT 



Aims. We model the chemistry of the inner wind of the carbon star IRC+10216 and consider the effects of periodic shocks induced 
by the stellar pulsation on the gas to follow the non-equilibrium chemistry in the shocked gas layers. We consider a very complete 
set of chemical families, including hydrocarbons and aromatics, hydrides, halogens, and phosphorous-bearing species. Our derived 
abundances are compared to those for the latest observational data from large surveys and the Herschel telescope. 
Methods. A semi-analytical formalism based on parameterised fluid equations is used to describe the gas density, velocity, and 
temperature from 1 Rj, to 5 R*. The chemistry is described using a chemical kinetic network of reactions and a set of stiff, ordinary, 
coupled differential equations is solved. 

Results. The shocks induce an active non-equilibrium chemistry in the dust formation zone of IRC+10216 where the collision 
destruction of CO in the post-shock gas triggers the formation of O-bearing species such as H2O and SiO. Most of the modelled 
molecular abundances agree very well with the latest values derived from Herschel data on IRC+10216. The hydrides form a family 
of abundant species that are expelled into the intermediate envelope. In particular, HF traps all the atomic fluorine in the dust formation 
00- zone. The halogens are also abundant and their chemistry is independent of the C/0 ratio of the star. Therefore, HCl and other CI- 



bearing species should also be present in the inner wind of O-rich AGB or supergiant stars. We identify a specific region ranging from 
2.5 K^, to 4 Rj,, where polycyclic aromatic hydrocarbons form and grow. The estimated carbon dust-to-gas mass ratio derived from the 
mass of aromatics formed ranges from 1.2 X 10"' to 5.8 x 10"' and agrees well with existing values deduced from observations. This 
Q ■ aromatic formation region is situated outside hot layers where SiC2 is produced as a bi-product of silicon carbide dust synthesis. The 

' MgS grains can form from the gas phase but in lower quantities than those necessary to reproduce the strength of the 30 //m emission 

, band. Finally, we predict that some molecular lines will show a flux variation with pulsation phase and time (e.g., H2O), while other 

d ■ species will not (e.g., CO). These variations merely reflect the non-equilibrium chemistry that destroys and reforms molecules over a 

pulsation period in the shocked gas of the dust formation zone. 

^SJ ■ Key words. Stars: Carbon - Astrochemistry - Stars: AGB and post-AGB - Molecular processes 

> 

] 1 . Introduction AI2O3). Conversely, if the star has experienced third-dredge up 

V ^ . in the upper part of its AGB ascension, its photosphere may be- 

VL; . In their late stages of evolution, low-mass stars (i.e., stars with ^^^^ carbon-rich. The chemistry of the wind then reflects the 

^ ■ initial masses on the Zero-Age-Main-Sequence comprised be- ^^^^^^ ^^^^^^ j^^j^^^ and is rich in C-bearing 

^ . tween 1 and 8 M^) ascend the Asymptotic Giant Branch (AGB) -^^ ^^^j^ acetylene, C2H2, and hydrogen cyanide, HCN. 

^ . and develop cool and strong stellar winds characterised by a ^j^^^^ ^^^^ j^g^ ^^^^^^^ ^^^^^^ ^^^^ ^j^^^ ^^^^-^ 

r-l ■ great variety of chemical species in the outflow detected through hotosphere. This simple picture was seriously questioned by 



X 



their ro-vibrational transitions (Olofsson 2008). With the launch detection at millimetre (mm) wavelengths of HCN in O-rich 
of the submillimetre (submm) Herschel telescope and the be- ^^^^ (Deguchi & Goldsmith 1985). In carbon stars, SiO 



ginning of science operations of the Atacama Large Millimetre observed (Olofsson et al. 1982), but because TE models of 

H . Array (ALMA), our knowledge of the chemical composition of ^^^^^^ ^^^^^ predicted its formation, it did not come as a surprise 

. ^ ■ AGB winds IS bound to dramatically increase with the discovery ^^^^^^^ ^^^-^ O-bearing species in carbon-rich environments, 

and Identification of many new molecules. The wind of a AGB However, the observed abundances were much higher than those 

star develops in the dense and hot gas layers above the stellar ^^^-^^^ ^^^^ S^^g^^j mechanisms were proposed to explain 

photosphere, triggered by the formation of dust. To a first ap- ^j^^^^ unexpected species, including an ion-molecule chemistry 

proximation, the chemical composition of the photosphere is de- ^j^^ ^^^^^ p^^^ ^j^^ ^^^jj^^ ^-^^ experiencing the penetra- 

termined by thermodynamic equilibrium (TE) owing to the high ^-^^ ^j^^ ultraviolet (UV) interstellar radiation field and cosmic 

temperatures and densities of the gas (Tsuji 1973, McCabe et ^^^^ ^^ejad & Miflar 1982). The later detection by the Infrared 

al. 1979). Carbon monoxide, CO, which is a very stable species, ^p^^^ Observatory, ISO, of hot vibrational transitions of carbon 

forms under TE conditions after the production of molecular hy- ^-^^-^^^ ^O., in the cool supergiant NML Cyg (Justtanont et al. 

drogen, H2. If the stellar photosphere is oxygen-rich, the excess ^995) indicated that carbon-bearing species could form in the 

oxygen not locked up in CO drives a wind chemistry dominated ^^^p j^^^^^ ^^^jj^^ ^-^^^ ^^^^ detection of water, H,0, 

by oxygen-bearing molecules such as water, H2O, and silicon ^-^^^ ^j^^ g^^S ^^^^jji^^ ^^^^^ wavelength by Melnick et'al. 

monoxide, SiO, and the dust forming in the wind includes sili- (2OOI) and that of hydroxyl, OH, by Ford et al. (2004) in the car- 
cates (e.g., forsterite Mg2Si04) and metal oxides (e.g., alumina 
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bon star IRC+10216 provided additional evidence of the com- 
plex chemistry of AGB outflows, where comets (Melnick et al. 
2001), grain-surface chemistry (Willacy 2004), or UV photodis- 
sociation in the intermediate envelope (Agiindez & Cernicharo 
2006) were proposed as possible sources of water. 

However, that these unexpected species could form by means 
of some non-equilibrium chemistry in the inner wind was soon 
found to be viable by surveys of specific species in a large sam- 
ple of objects, e.g., SiO (Schoier et al. 2006). Just above the pho- 
tosphere, the gas experiences the passage of shocks triggered by 
the pulsation of the star and these shocked gas layers are the 
locus of dust formation (Bowen 1988). These dense molecu- 
lar gas layers has been detected by both observations with ISO 
(Tsuji et al. 1997) and near-infrared (IR) interferometry (Perrin 
et al. 2004). A semi-analytical model for the physics of these 
shocked regions based on the work of Fox & Wood (1985) and 
Bertschinger & Chevalier (1985) was proposed by Cherchneff 
(1996) assuming that the shock energy was dissipated in the im- 
mediate postshock gas by the collisional dissociation of H2 and 
that subsequent cooling was provided by adiabatic expansion. 
Using a similar formalism for the inner wind of IRC+10216, 
Willacy & Cherchneff (1998) modelled the gas-phase chem- 
istry and showed that the collisional destruction of CO in the 
shocks could release free atomic oxygen and trigger the forma- 
tion of OH and SiO. The formations of HCN, CO2, and CS were 
later described in the inner wind of the O-rich Mira star IK Tau 
by Duari et al. (1999). These two studies therefore highlighted 
the importance of the shock-induced non-equilibirum chemistry 
in the inner wind of AGB stars in unleashing the synthesis of 
molecules that were not expected to form under photospheric 
TE conditions. More generally, Cherchneff (2006) studied the 
inner wind composition as a function of the carbon-to-oxygen 
(C/O) ratio of the stellar photosphere and confirmed the forma- 
tion of O-bearing species in carbon stars and that of C-bearing 
molecules in O-rich AGB stars as a result of non-equilibrium 
chemistry induced by periodic shocks. Several observations of 
high energy rotational transitions of HCN, SiS, CS, and SiO in 
AGBs and supergiants at mm and submm wavelengths corrobo- 
rated these results (e.g., Schoier et al. 2006, 2007, Ziurys et al. 
2007, 2009, Decin et al. 2008). 

A next step in understanding the complexity of the wind 
chemistry close to the stellar photosphere was achieved with 
Herschel and the confirmation of the widespread presence of wa- 
ter in the dust formation zone of several carbon stars (Decin et 
al. 2010b, Neufeld et al. 2010, 2011a, 2011b). Agundez et al. 
(2011) proposed that the partial dissociation of '^CO and SiO 
by UV photons penetrating deep inside the clumpy stellar wind 
could produce free atomic oxygen and lead to the subsequent 
formation of H2O. This scenario would destroy other species 
in the deep layers and have some impact on the water isotopo- 
logue abundances as discussed by Neufeld et al. (2011a). Both 
effects have not yet been tested observationally. Updating the 
inner wind model for IRC+10216, Cherchneff (2011a) showed 
that the non-equilibrium chemistry triggered by the shocks could 
form water in the dust formation zone in competition with the 
synthesis of SiO, with abundances that agree excellently with 
those derived from observations. 

In this study, we present a complete, updated chemical model 
of the inner wind of the carbon star IRC+10216, based on the 
non-equilibrium chemistry approach mentioned previously. The 
synthesis of classical molecules including water and several new 
species is considered in the gas-phase, comprising hydrides, 
chlorine and phosphorous compounds, hydrocarbons and aro- 
matics, and gas-phase precursors to silicon carbide and metal- 



Table 1. Stellar parameters for IRC+10216 



(Mo) 


2 


1 


R*(cm) 


6.5xl0'3 


2 


C/O 


1.4 


3 


P (days) 


650 


4 


r. (R*) 


1.2 


2 


V, (km s-') 


20 


2 




3.6xl0'3 


5 


T,„.(r,) (K) 


2062 


5 



References. (1) Willacy & Cherchneff 1998; (2) Ridgway & Keady 
1988; (3) Winters et al. 1994; (4) Witteborn et al. 1980; (5) Cherchneff 
2011a. 



sulphide dust. The goal of such a study is to confirm the for- 
mation of already detected species and predict new molecules 
potentially observable by means of their high excitation transi- 
tions in the inner wind. In Section|2l we give a brief description 
of the physics of the shocked layers, while the chemistry is dis- 
cussed in Section[3] the results for the various families of species 
are presented in SectionH] and a discussion follows in Section|5] 

2. Physics of the inner wind 

The stellar parameters used in this study for IRC+10216 are 
listed in Table [T] The stellar photosphere and the outer atmo- 
spheric layers above it are assumed at TE with a solar elemental 
composition (Asplund et al. 2009) except for carbon where the 
C/O ratio was set to 1.4 (Winters et al. 1994). The shocks are 
assumed to form at a radius rj = 1.2 R* with a velocity V,. The 
TE calculations are run until rj, which is characterised by the 
temperature and number density listed in Table [T] and the TE 
abundances are used to characterise the unshocked gas. Since the 
pre-shock gas at r, has a high molecular component, the impact 
of the passage of periodic shocks is modelled by considering a 
post-shock gas that cools via both 1) the collision-induced dis- 
sociation of H2 and 2) adiabatic expansion. The shock jump in 
density, temperature, and velocity is described by the Rankine- 
Hugoniot jump conditions applied to the pre-shock gas at r^ for 
a shock velocity = 20 km s"' (Ridgway & Keady 1981). 
According to Fox & Wood (1985), radiative processes do not 
operate in cooling the post-shock gas owing to the modest shock 
strength, for which the energy loss is provided by the endother- 
mic dissociation of H2 by collisions initiated by the prevalent 
chemical reaction 

H2 +H2 ^ H + H + H2. (1) 

Reaction [T] has a reaction rate given by kdiss - 1.5 x lO"** x 
exp(-48346/T) cm^ s"', where T is the gas temperature, and 
operates over a length 1^,^^ defined as 



1 



kdiss X n(H2) 



N;, 



(2) 



where n(H2) is the number density of H2 in the post-shock gas, 
kdiss is the rate of reaction [1] V , is the shock velocity at r, as 
given in Table [1] and Ny,,,,,,, is the Rankine-Hugoniot velocity 
shockjump. 

Once this cooling has proceeded over 1^,5^, we assume that 
later cooling occurs via adiabatic expansion. As modelled by 
Bertschinger & Chevalier (1985), the fluid equations that de- 
scribe the conservation of mass, momentum, and energy are 
parametrised and solved for the boundary conditions imposed 
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by stellar gravity and the return of the shocked gas to its ini- 
tial pre-shock position. Typical excursions of the gas layers over 
several pulsation periods are illustrated in Willacy & Cherchneff 
(1998) and Cherchneff (2011a) for IRC+10216. The pre-shock 
gas temperature and density profiles are derived as a function of 
radius using the formalism of Cherchneff et al. (1992), where the 
impact of shocks on the gas density is described by an extended 
scale-height formalism for the initial conditions listed in Table 
[1] The derived pre- and post-shock gas parameters at various ra- 
dius are those assumed in Cherchneff (201 la) and listed in Table 

12 

3. Chemistry of the inner wind 

New families of chemical species are considered along with 
the classical molecules already studied in previous models 
(Cherchneff 2006, 2011a). They include hydrides, halogens, 
in particular chlorine- and fluorine-bearing species, and finally 
phosphorous-bearing species. Small gas-phase molecular pre- 
cursors to silicon carbide and metal-sulphide grains are also 
included. All species are listed in Table [3] where the largest 
molecule is the aromatic ring of benzene, CeHe. For metal hy- 
drides, there is no available information on chemical reaction 
rates for the formation processes of several of them, namely 
FeH, MgH, NaH, KH, and PH, at the high temperatures char- 
acterising the inner wind. We thus assume rates similar to those 
of documented reactions using the principle of isovalence as a 
guideline. Generally, the typical rates of the dominant formation 
pathways that consists of the reaction of the metal with H2 have 
an Arrhenius factor ranging from 10"'^ and lO"** cm^'s"', an ac- 
tivation energy barrier of a few 1000 K, and a mild temperature 
dependance (e..g, Cohen & Westberg 1983). 

The chemistry of halogens including chlorine- and fluorine- 
bearing species is rather well-documented at high and intermedi- 
ate temperatures. However, this is not the case for the phospho- 
rous chemistry. Owing to the isovalence of phosphorous, P, with 
nitrogen, N, we assume for the P-bearing species chemical pro- 
cesses similar to those involving N and for which the rates are 
measured or calculated. Finally, the chemistry of silicon. Si, and 
sulphur, S, is poorly known and studied. For these elements, we 
restrict the chemical processes to the set of reactions for which 
reaction rates are documented. 

All chemical pathways that lead to the formation of linear 
molecules, carbon chains, and aromatic rings, include neutral- 
neutral processes such as termolecular, bimolecular, and radia- 
tive association reactions, whereas destruction is described by 
thermal fragmentation and neutral-neutral processes (i.e., oxi- 
dation reactions of hydrocarbons and all reverse processes of 
the formation reactions). No ions are considered in this chem- 
istry because the UV stellar radiation field of IRC+10216 is too 
weak to foster the efficient photodissociation and ionisation of 
molecules. In total, 63 species were considered for a chemical 
network of 377 reactions. Details of all these processes are pro- 
vided in the online Appendix 1 available at the CDS, which gives 
the reactions included in the chemical network and their reaction 
rates. It contains the following information: Column 1 lists the 
reaction number. Column 2 gives the reactants and the products. 
Column 3 lists the A coefficient. Column 4 lists the n factor. 
Column 5 lists the activation energy in Kelvin, and Column 6 
gives the reference for the data. The major differences between 
the present chemical network and that used in previous stud- 
ies (e.g., Willacy & Cherchneff 1998, Cherchneff 2006) are ex- 
plained in detail in the Appendix of Cherchneff (201 la). In terms 
of formalism, three major changes are implemented. Firstly, the 



treatment of the reverse reaction of a specific process is changed. 
Several new rates have been measured in combustion and aerosol 
chemistry, and are now available. The calculation of the reverse 
rate from the equilibrium constant assuming detailed balance 
(i.e.. Equation 4 in Willacy & Cherchneff 1998) can lead to er- 
roneous values when the gas temperature and density decrease. 
Therefore, we prefer to directly enter the available measured or 
calculated rate values in the chemical network. When the infor- 
mation is unavailable, we make 'educated' guesses depending 
on the type of the reaction and its endo-(exo)thermicity. This ap- 
proach allows us to test a variety of reaction paths that would 
be closed if one were to strictly apply detailed balance consid- 
erations. Secondly, the chemistry involving atomic silicon. Si, 
and Si-bearing species is restricted to reactions for which rates 
have been measured or calculated. Thirdly, as previously men- 
tioned, the chemistry now describes the formation of a larger set 
of chemical species. This chemical scheme is used to solve a 
set of 63 stiff, coupled, ordinary differential equations (ODEs) 
at each radius of the inner wind. These coupled ODES are inte- 
grated over space (for the H2 cooling region defined by Equation 
|2]l and time (for the adiabatic expansion over a pulsation period) 
for the radii and corresponding shock strengths of Table |2] The 
post-shock abundances of the species at the end of the pulsation 
period and radius r, are used as pre-shock initial abundances for 
the successive radius r,+i. 

In all, the chemical network aims to provide a comprehensive 
and coherent chemical description of the wind, where known 
trends are reproduced or new trends are presented. It does not 
aim to give an exact description of all chemical processes in the 
wind, which is an impossible task to perform. Owing to the un- 
certainties in both our theoretical description of the inner wind 
and observational data analysis, we consider in the following 
sections that a good agreement has been reached when modelled 
abundances and derived values from observations differ by at 
most a factor of ten. 



4. Results 

Results for the dominant chemical species in the inner wind are 
summarised in Table|4] We list the abundance values derived by 
applying TE to the gas conditions met at r,, as well as the non- 
equilibrium abundances at r, and 5 R*, and abundances derived 
from the most recent observations. As stated before, all mod- 
elled and observational values differ by at most a factor of ten, 
except for PN (see 34.41 ). This good agreement is illustrated in 
Figure [1] where we compare the modelled and observed abun- 
dances with respect to H2 at 5 R*. Discrepancies between values 
derived from TE and non-equilibirum chemistry for O-bearing 
species are clear from inspection of Table |4] These discrepan- 
cies highlight once more the importance of the shock chemistry 
in the dust formation zone of AGB stars. We discuss below the 
results for the specific chemical families under study. 

4. 1 . Prevalent molecules 

A group of species dominates the molecular phase of the 
shocked, inner wind along with H2. It consists of CO, HCN, 
CS, N2, C2H2, SiS, SiO, and H2O, whose abundances relative 
to H2 are illustrated as a function of radius in Figure |2] while 
their abundances with respect to total gas and as a function of 
pulsation phase 9 at the shock formation radius r , are shown in 
Figure[3] All molecules experience destruction in the post-shock 
gas at 6 < 0.05 as seen in Figure [3] The destruction is more 
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Table 2. Gas parameters in the dust formation zone of IRC+10216 taken from Cherchneff (2011a). The shocks form at r, = 1.2 R*. 
For each radius, the temperature and number density are given for the pre-shock gas, the gas after the shock front (in the collision- 
induced H2 dissociation zone), and the gas at the beginning and the end of the adiabatic cooUng zone. The temperature T is in Kelvin 
and the gas number density ngas in cm"^ . 



r 


Shock strength 


Preshock gas 


Shock Front 


Adiabatic 


expansion - start 


Adiabatic 


expansion - end 


(RJ 


(kms-') 


T 


^gas 


T 


^gas 


T 




T 




1.2 


20.0 


2062 


3.63(13) 


19725 


1.98(14) 


4409 


5.97(14) 


1480 


3.63(13) 


1.5 


17.9 


1803 


8.24(12) 


15922 


4.40(13) 


3870 


1.29(14) 


1290 


8.24(12) 


2 


15.5 


1517 


1.44(12) 


12081 


7.59(12) 


3200 


2.14(13) 


1080 


1.44(12) 


2.5 


13.9 


1327 


4.24(11) 


9779 


2.21(12) 


2750 


6.08(12) 


951 


4.24(11) 


3 


12.6 


1190 


1.69(11) 


8245 


8.73(11) 


2430 


2.35(12) 


848 


1.69(11) 


4 


11.0 


1001 


4.51(10) 


6284 


2.29(11) 


1790 


4.48(11) 


711 


4.51(10) 


5 


9.8 


876 


1.79(10) 


5096 


8.94(10) 


1550 


1.71(11) 


621 


1.79(10) 



Table 3. Chemical species included in the inner wind model of IRC+10216. 
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Fig. 1. Observed and modelled abundances (with respect to H2) 
for most of the molecules listed in Table |4] Abundances are at 
phase - I and 5 R*. The agreement between modelled and 
observed values is satisfactory and less than a factor of ten for 
all species apart from PN (for CO both values overlap). 



Fig. 2. Abundances (with respect to H2) of the prevalent 
molecules formed in the inner wind as a function of radius. 
Abundances are those at phase 6=1, i.e., after one full oscil- 
lation cycle. Abundances derived by assuming thermodynamic 
equilibrium at r^ are moved to radius 1 R* for clarity. 



or less severe depending on the species and no large discrepan- 
cies may exist between TE and non-TE abundances for specific 
molecules at the end of one oscillation. However, TE abundance 
values for several species (CS, SiS, SiO, and H2O) differ drasti- 
cally from those obtained from shock-induced, non-equilibrium 



chemistry, as already stressed in existing studies (Willacy & 
Cherchneff 1998, Cherchneff 2006, 201 la). 

The formation of O-bearing species, namely H2O and SiO, 
results from the collision dissociation of CO in the post-shock 
gas. For the specific conditions of the IRC+10216 model, be- 
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Table 4. Predicted abundances (with respect to H2) of selected species. TE abundances are given at r, and modelled abundances are 
listed for r, and 5 R*. Most recent abundances (with respect to H2) derived from observations are also listed when available with 
the corresponding reference. W corresponds to abundances lower than 1 x 10 



Species 


TE at r, 


Predicted at r, 


Predicted at 5 R^, 


Observed 


Reference 


CO 


1.5 X 10-^ 


8.1 X 10- 


-4 


9.1 X 10-" 


1.0 X 10-^ 


Herschel - Cernicharo et al. (2010a) 


HCN 


4.3 X 10-5 


3.2 X 10- 


-5 


5.3 X 10-'^ 


8.0 X 10-^ - 2.0 X 10-5 


IRAM - Cernicharo et al. (201 1) 


C2H2 


9.2 X 10-5 


2.0 X 10- 


-4 


1.5 X lO-" 


8.5 X 10-5 


TEXES/IRTF - Fonfn'a et al. (2008) 


N2 


7.5 X 10-5 


4.4 X 10- 


-5 


5.8 X 10-5 






SiS 


4.3 X 10-« 


2.2 X 10- 


-5 


2.1 X 10-5 


2.0 X 10-« - 8.0 X 10-* 


PACS/SPIRE Herschel - Decin et al. (2010a) 


CS 


3.6 X 10-5 


5.6 X 10 


-6 


6.4 X 10-'^ 


9.3 X 10-*^ 


SMA - Patel et al. (2009) 


SiO 


4.1 X 10-' 


8.4 X 10- 


-6 


3.6 X 10-** 


2.0 X 10-** -3.0 X 10-' 


PACS/SPIRE Herschel - Decin et al. (2010a) 


H2O 


2.0 X 10-'- 


9.9 X 10 


-5 


7.1 X 10-' 


8.0 X 10-** 


HIFI/Herschel - Neufeld et al. (201 la) 




- 


_ 






> 1 X 10-'° 


PACS/SPIRE Herschel - Decin et al. (2010b) 


SiCj 


4 X 10-^ 


3.8 X 10- 


11 


5.5 X 10-** 


2.0 X 10-' 


HIFI/Herschel - Cernicharo et al. (2010a) 


HCl 


5.4 X 10-^ 


6.3 X 10 


-7 


4.9 X 10-' 


1.0 X 10-' 


HIFI/Herschel - Agundez et al. (201 1) 


HF 


1.2 X 10-'' 


7.3 X 10- 


-8 


7.2 X 10-** 


8.0 X 10-' 


HIFI/Herschel - Agundez et al. (201 1) 


NaCl 


2.4 X 10-'^ 


6.7 X 10- 


-9 


1.2 X 10-** 


2 X 10-' 


SMT/ARO - Milam et al. (2007) 


AlCl 


1.6 X 10-'" 


1.8 X 10- 


10 


1.3 X 10-' 


3.5 X 10-** 


IRAM - Agundez (2009) 


KCl 


6.4 X lO-''' 


7.2 X 10- 


10 


1.2 X 10-' 


7.0 X 10-'" 


IRAM - Agundez (2009) 


AlH 


3.8 X 10-' 


4.6 X 10 


-6 


2.0 X 10-** 


< 1.6 X lO-*-" 


PACS/SPIRE Herschel - Cernicharo et al. (2010b) 


NaH 


1.1 X 10-" 


1.4 X 10- 


-6 


9.9 X 10-'2 


< 3.0 X lO-''' 


PACS/SPIRE Herschel - Cernicharo et al. (2010b) 


KH 


2.4 X 10-'- 


1.1 X 10 


-7 


9.0 X 10-'^ 


< 3.0 X 10"'" 


PACS/SPIRE Herschel - Cernicharo et al. (2010b) 


MgH 


3.7 X 10-' 


7.9 X 10- 


-9 
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Fig. 3. Abundances (with respect to total gas) of the prevalent 
molecules shown in Figure|2]as a function of phase 9 at the shock 
formation radius r,. Destruction occurs at small values of 9 but 
the molecules reform as the gas cools down. 



tween 10 % and 20 % of CO molecules are destroyed at r, in 
the H2 dissociation cooling part, while they quickly reform in 
the adiabatic expansion at phase 9 > 0.05, making CO the main 



provider of atomic oxygen in the post-shock gas. The formation 
of H2O is thus triggered by the reactions (Cherchneff 201 la) 



O -H H2 ^ OH -H H 
and 

OH + H2 ^ H2O + H 



(3) 



(4) 



The reaction given in Eq.|4]is in competition with the formation 
of SiO via the reaction 



Si -H OH ^ SiO + H. 



(5) 



The rate of the backward process of the reaction in Eq.|5]is un- 
known but the reaction has an endothermicity of ~ 40000 K at 
2000K. As mentioned in Section[3] we assume a low rate for that 



reaction (k = 1 x 10 



s ) that reflects its low efficiency 



over the temperature range of interest. We considered the result- 
ing water abundance for different rate values and temperature 
dependences. The water abundance always shows a trend sim- 
ilar to that reported here, i.e, a high inner value that decreases 
at radius r > 2 R* to reach a typical value of ~ 1 x lO"^. In 
contrast, the rate of the backward reaction in Eq. |5]was explic- 
itly calculated from the equilibrium constant given in the study 
of Willacy & Cherchneff (1998). This previous rate had a very 
low value for temperatures lower than 2000K, and contributed 
in part to the non-replenishment of OH and the disappearance of 
H2O at larger stellar radii. 



5 



I. Cherchneff: The dust formation zone of IRC+10216 revisited 



1x10"''° 



1 

- IN 

>-^'. — ^'•r: — 


AIH 


HCI 




I; 


; ^ 

HF 






\ 


/V '""••••...NaH 






1 if 

1 1SH\ 


\ MgH ^\ 
FeH\ ^ \, 

1 \. 1 








2 3 

Radius (Rt) 


4 





Fig. 4. Abundances (with respect to H2) of the hydride species 
formed in the inner wind as a function of radius. Abundances 
derived by assuming thermodynamic equilibrium at r., are moved 
to radius 1 R* for clarity. 

The prevalent molecules in Figure |2] come from different 
chemical families that are chemically linked together by the indi- 
rect key-role of the over-abundant H2 species. Molecular hydro- 
gen is primarily involved in the formation processes of both vari- 
ous members of chemical families (e.g., hydrides, hydrocarbons) 
and specific molecules such as OH and H2O. The destruction or 
formation of H2 can then impact all chemical families, which be- 
come interrelated. This is particularly true for water, owing to its 
link to the hydrocarbon family. The destruction of hydrocarbons, 
starting with C2H2, indeed releases H2. Through the synthesis of 
hydroxyl, OH, via the reaction of atomic O with H2 in the post- 
shock gas, H2O is thus linked to the hydrocarbon family. As seen 
above, water is also linked to the Si-bearing species as both H2O 
and SiO species are competitors in the depletion of OH. Atomic 
Si is also efficiently included in silicon monosulphide, SiS, and 
in silicon carbide, SiC, in the hot post-shock gas at r,. Therefore, 
the Si-bearing species through their connection to SiO impact 
the water abundance, whose value listed in Table |4] is slightly 
higher than the value derived by Cherchneff (201 la). The chem- 
istry of SiC dust precursors has been extended to the first ring 
(SiC)2 in the present study, and as discussed in § 14.61 the forma- 
tion of SiC and the rings SiC2 and (SiC)2 proceeds very early on. 
The more atomic Si is trapped in SiC, the less Si is able to react 
with OH to form SiO. Thus, the SiO abundance listed in Table 
|4]is lower than that derived by Cherchneff (2011a) by a factor 
of about two, resulting in a higher water abundance. However, 
the chemical trends and processes in both studies are alike. The 
formation pathways to other important molecules such as C2H2 , 
CS, and HCN were discussed in detail by Cherchneff (2006) for 
a carbon star with a C/O ratio equal to 1 . 1 and similar chemical 
routes operate in the inner wind of IRC+10216. Overall and as 
seen from Figure [1] the most abundant species have modelled 
abundances that agree very well with values derived from obser- 
vations. 

4.2. Hydrides 

The interest of studying light hydrides was rekindled with the 
launch of Herschel. The rotational spectra of these species lie in 
the submm and far-IR domains and are difficult to observe from 
Earth. The detection by HIFI onboard Herschel of the J=l-^0, 



2— >1 and 3— »2 transitions of HCI and the J=l— >0 transition of 
HF in IRC+10216 was reported by Agundez et al. (201 1). 

The gas-phase chemistry of many light hydrides has never 
been clearly characterised and only a few measured rates have 
been documented. However, their production in the laboratory 
from the gas phase occurs via the reaction of a metal vapour with 
hydrogen (Ozin & McCaffrey 1984). For our model, we con- 
sider all known chemical reactions and extend the documented 
chemistry to some species according to the isovalence of specific 
elements (e.g., N and P). Following the prescriptions of exper- 
imental studies, we model the formation of hydrides according 
to 

X + H2 ^ XH + H, (6) 

where X represents any atomic species. H2 is destroyed in the 
strongest shocks up to ~ 3 R* but efficiently reforms in the post- 
shock gas at phase 6 > 0.3 (see Figure[3]l. The large H2 reservoir 
in the gas layers insures that most hydrides are formed follow- 
ing the reaction in Eq. |6] The abundance variation with radius 
of the major hydrides is illustrated in Figure H) Hydrides have 
variations that depend on the species. The abundances of AlH 
are high out to 3 R* but decrease at larger radii owing to the for- 
mation of aluminium chlorine, AlCl (see 34.3l below). A similar 
behaviour applies to NaH and KH with the formation of NaCl 
and KCl. In contrast, HCI and HF have consistently high abun- 
dances (a-(HCI) ~ 4 X 10"^ and x(HF) ~ 7 x 10"**) in the inner 
wind while MgH, FeH, PH, and SH have extremely low abun- 
dances that are well below the PACS/SPIRE 3cr detection limits 
estimated by Cernicharo et al. (2010b). 

The Herschel detection of the J=l— >0 transition of HF by 
Agundez et al. (2011) points to a constant abundance of ~ 
8 X 10"^ with respect to H2 extending from the inner envelope up 
to 45 R* . This value is lower by a factor of nine than that derived 
from TE in the photosphere, where most of the fluorine is in the 
form of HF. To reconcile these two values, they claim that F must 
be depleted on dust grains in the inner wind. Our calculated HF 
abundance in Table|4]has the constant value of 7.25 x 10"^, cor- 
responding to the solar abundance for F. This value reflects the 
quick conversion of fluorine into HF by its reaction with H2 in 
the post-shock gas at all radii in the dust formation zone ex- 
tending to 5 R*. The result confirms that HF acts as the main 
reservoir of fluorine in AGB stars. The discrepancy between our 
kinetic results and those of Agundez et al. is difficult to quan- 
tify. Our abundance variation through the inner wind agrees well 
with the abundance profile that is required to reproduce the HF 
HIFI data, i.e., a constant abundance distribution extending up to 
45 R*, but our calculated abundance value is higher Agundez et 
al. quote that the error in the radiative transfer model is a factor 
of two, while the reaction rate for the gas-phase formation of HF 
from H2 has been measured and is well-documented. The use 
of other available rate values does not change the present result, 
which indicates that fluorine is totally depleted into hydrogen 
fluoride in the dust formation zone. 

4.3. Chlorines 

Chlorine (CI) has a solar abundance with respect to hydrogen 
of ~ 3.2 X 10"^ (Asplund et al. 2009), and is present at TE in 
the photosphere in the form of both HCI and CI of comparable 
abundances with respect to H2 of ~ 2.6 xlO"^. Chlorine-bearing 
species have long been observed in circumstellar envelopes of 
AGB stars. Recent large molecular surveys confirm the presence 
of NaCl, AlCl, and KCL in IRC+ 1 02 1 6 (Tenenbaum et al. 20 1 0) 
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Fig. 5. Abundances (with respect to H2) of chlorine- and 
phosphorus-bearing species as a function of radius. Abundances 
derived by assuming thermodynamic equilibrium at rj are moved 
to radius 1 R* for clarity. 

when the hydride HCI was detected with Herschel and unam- 
biguously found very close to the star (see § 14. 2b . Shinnaga et al. 
(2009) also identified KCI lines in their e-SMA observations that 
have a very compact distribution centred on the star Results for 
the prevalent CI -bearing species are shown in Figure |5] Clearly, 
HCI is still the dominant Cl-bearing species in the dust forma- 
tion zone, and readily forms at rj with an almost constant abun- 
dance throughout the inner wind and a value of 3.9 x 10"^ at 5 
R* . This values perfectly agrees with that derived by Agundez 
et al. (201 1) and is very close to the TE abundance of HCI at r,. 
However, HCI is destroyed in the shock at rj and reforms in the 
post-shock region at phase Q > 0.2 to finally reach an abundance 
after one oscillation very close to the TE value. Other Cl-bearing 
species, namely AlCl, NaCl, and KCI in decreasing order of im- 
portance, are also present in the wind and their production orig- 
inates directly from that of HCI via the reaction 

X H- HCI ^ XCl + H, (7) 

where X is the metal. All documented reactions of Eq. Q are 
fast, with activation energy barriers of a few thousand degrees 
(e.g., Husain & Marshall 1986), requiring high temperatures to 
proceed. 

Very much like H2 for hydrides, HCI acts as the production 
agent of Cl-bearing species in the dust formation zone of AGB 
stars. The chemistry of chlorine depends essentially on the hy- 
drogen and the chlorine content of the gas through HCI, and is 
thus independent of the C/O ratio of the photosphere. Since HCI 
is a rather stable molecule (with a dissociation energy Do - 4.4 
eV), it should also form efficiently in oxygen-rich AGB stars and 
by a chemical pathway similar to the reaction in Eq. Rowing to 
the large H2 reservoir of AGB winds. Other Cl-bearing species 
will also form via the reaction in Eq.|2]provided that metal atoms 
are available in the gas phase. HCI and other Cl-bearing species 
are thus expected to be present with intermediate-to-high abun- 
dances in the inner wind of oxygen-rich AGB stars. The detec- 
tion of high energy transitions of NaCl towards the O-rich su- 
pergiant VY CMa and the O-rich Mira IK Tau by Milam et al. 
(2007) indeed indicates that these two species form close to the 
star, with abundances of 5x10"^ and 4x10"'*, respectively. Our 
results corroborate these observations and we predict that HCI, 
NaCl, and KCI should be observable in the dust formation zone 



of O-rich evolved stars. The abundance of AlCl may be lower 
because 1) the Al -bearing species AlOH has a high abundance 
in the wind acceleration zone - Tenenbaum & Ziurys (2010) de- 
rived a AlOH abundance of ~ 1 x 10"^ for VY Cma - and 2) a 
large fraction of Al is expected to be depleted in alumina, AI2O3, 
in the stellar wind. In reality, AlCl in VY CMa was not detected 
by Tenenbaum & Ziurys (2010). 

4.4. Phosphorous bearing species 

The phosphorous-bearing molecules HCP, CP, C2P, and PN have 
been detected in the wind of IRC+10216 (Guelin et al. 1990, 
Agundez et al. 2007, Milam et al. 2008, Halfen et al. 2008, He 
et al. (2008), Tenenbaum et al. 2010). While the shapes of the 
C2P and CP line profiles are indicative of a shell-like distribution 
with a formation locus in the outer envelope induced by U V pho- 
todissociation, HCP and PN have been found close to the star. 
The chemistry of phosphorous is poorly documented and we use 
the isovalence of P with N to estimate the rates of a set of basic 
formation and destruction processes derived from the equivalent 
processes involving N. We consider a few molecules including 
PN, HCP, CP, and P2. The latter species was chosen to reflect the 
refractory nature of phosphorous and its ability to form clusters. 
Results for P-bearing molecules are shown in Figure|5] PN is the 
prevalent phosphorous compound followed by HCP and P2. PN 
is chiefly formed at r s by the two reactions 

N + CP ^ PN + C (8) 

and 

CN + CP ^ PN + C2 (9) 

Similar reactions for the isovalent element N have measured 
rates at high temperatures. The resulting abundances for PN are 
quite high (~ 3 x 10"^) through the inner wind with a rapid for- 
mation at Ts where the molecule reaches its final abundance at 
phase = 0.5 in the post-shock excursion. This is coherent with 
the fact that PN mimics N2 which shows a similar behaviour: a 
rapid synthesis in the post-shock gas at r, and a constant high 
value across the inner wind (see Figure|2]). At r, and at the early 
phases of the post-shock adiabatic excursion, phosphorous is 
quickly integrated into CP which later distributes P into HCP 
and PN. Milam et al. (2008) observed several rotational lines of 
PN and HCP in the inner wind of IRC+10216 and derived a low 
PN abundance of ~ 3 x 10 '". Our modelled value is greater 
by a factor of 1000. However, our derived HCP abundance of 
8.2 X 10"^ agrees well with the observations by Agundez et al. 
(2007) and Milam et al. (2008). The former study derived an 
abundance value of 1.4 x 10 for radii larger than 20 R*, and 
claimed that a depletion onto dust grains is necessary to recon- 
cile the high abundances derived assuming TE in the dust forma- 
tion region with those at 20 R*. The present results indicate that 
HCP depletion in dust may not be necessary because the high 
TE abundances at r < r, (~ 3.0 x 10"^) quickly drop to 8.2 x 10""* 
at larger radii owing to the non-equilibrium chemistry and the 
partial conversion of P into CP and PN in the post-shock gas. 

The high abundances of PN obtained in the model compared 
to observations force us to test the P chemistry and assess the 
conditions for which the conversion of CP into HCP and PN is 
effective. We decrease all chemical rates by a factor of ten but 
have no success in diminishing the PN abundance in the inner 
wind. We also lower the rates of the two forward reactions given 
in Eqs. |8]and|9] by a factor of 100, which results in decreasing 
the PN abundance to ~ 3.5 x 10"^ but in increasing the HCP 
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abundances to 3.2 x 10"^, which is high to agree with observa- 
tions. The assumption of isovalence between P and N requires a 
phosphorous chemistry in which HCP and PN mimic HCN and 
N2. In the inner wind, nitrogen is distributed between these two 
species in almost equal amounts, a result corroborated by the ex- 
cellent agreement of HCN abundances with values derived from 
observations. In the case of P, this distribution is not observed 
and the low PN abundance derived by Milam et al. (2008) indi- 
cates a low efficiency channel for the conversion of CP into PN 
in the dust formation zone. According to this result, they deduce 
that if PN/N2 were approximately equal to P/N, the Nt abun- 
dance would be very low (~ 1 x 10"^). Such a low value con- 
tradicts the high N2 abundance given by TE in the photosphere 
and fostered by the non-equilibrium chemistry of the inner wind 
(see Table m. The discrepancy regarding the abundances of PN 
questions the validity of the assumption of isovalence in the case 
of phosphorous, highlights the different chemical processes that 
may control the phosphorous chemical family, and points to the 
need for more high resolution observations of high energy tran- 
sitions of PN in the inner wind. 

Finally, inspection of Figure|5]reveals the presence of a mod- 
est amount of phosphorous dimers, P2, synthesized in the inner 
wind with x(P2) ~ 3 x 10"^ at 5 R*. As P has even lower abun- 
dances than P2, a small population of P clusters may grow from 
P2 coalescence but the aggregation process will terminate with 
the formation of the stable tetrahedral P4 cluster, as observed in 
laser ablation of red phosphorous crystals (Bulgakov al. 2002). 
At most, a P4 abundance of ~ 1x10"'" may form with a left-over 
population of dimers with abundances ~ 1 x 10 Therefore, 
phosphorous clusters should not be prevalent condensates in the 
dust formation zone of carbon stars. 



4.5. Carbon dust precursors: Hydrocarbons and aromatics 

The formation of the first aromatic ring of benzene represents 
a bottleneck to the formation of polycyclic aromatic hydrocar- 
bon (PAH) species and their growth. In the chemical scheme, 
it is described by the recombination of two propargyl radicals, 
C3H3 which is the dominant closure pathway, and the reaction 
of buten-3-ynyl radicals, C4H3, with C2H2. These two routes are 
the prevalent channels to aromatic formation in sooting flames 
on Earth (Cherchneff 2011b). The formation of C3H3 results 
from the reaction of C2H2 with methylene, CH2, in the imme- 
diate post-shock region at r,. Abundances with respect to H2 are 
shown in Figure |6] for the 20 km s ' shock at r,. At the high 
post-shock gas temperatures, only stable hydrocarbons such as 
C2H2 , C3H2, and C3H3 can form in large amounts. In Figure|7] 
the abundances for similar species are shown for a shock strength 
of 12.6 km s"' at 3 R*, with the appearance of C(,H(, at phases 

> 0.4. As apparent in Figure |2] once C2H2 forms at r^, it stays 
abundant over the inner wind region, providing a large reservoir 
to grow hydrocarbons, specifically C3H2. Once C3H3 starts to 
form from the reaction of C3H2 with H2, CeHe quickly builds 
up when the gas temperature drops. Indeed, higher temperature 
values favour the reverse routes to the formation of CeHs, which 
are endothermic channels. The clincher to build up the CgHe ring 
are thus the lower gas temperatures encountered at phases > 0.4. 
The abundances of hydrocarbons and aromatics as a function 
of radius are illustrated in Figure [8] The formation of aromat- 
ics is delayed to r > 2.5 R* because oxygen-bearing species are 
present at smaller radii. Water and hydroxyl are the main ox- 
idation agents to hydrocarbons and aromatics in the gas, and 
CfiHe forms when the O-bearing species abundances drop at r 

> 2.5 R*(see Figure |2]i. According to Figure [8] the prevalent 




Fig. 6. Abundances (with respect to H2) of hydrocarbon species 
as a function of pulsation phase 6 for the 20 km s"' at rj = 1.2 
R*. No aromatic forms at these small radii owing to the presence 
of the oxidation agents H2O and OH. 
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Fig. 7. Abundances (with respect to H2) of hydrocarbon and aro- 
matic species as a function of pulsation phase for the 12.6 km 
s ' at 3 R*. Benzene starts forming at the low gas temperatures 
(~ 1000 K) encountered at 6* > 0.4. 



hydrocarbon species that escape the inner envelope are C4H4, 
C4H2, and C3H3, although the growth and condensation of aro- 
matics to amorphous carbon (AC) grains may alter this result to 
some extent. 

Central to this study is the formation of the first aromatic 
ring as a bottleneck process to AC dust synthesis in carbon stars. 
As seen in Figure [8] the reported abundance of CeHg gradually 
increases with stellar radius up to ~ 3 R*, where it reaches a 
maximum value to gradually decrease up to 5 R*. This result 
points to a specific radius range where AC dust has a chance 
to efficiently grow provided that the growing agent, in this case 
acetylene, C2H2, is still abundant in the gas phase. To test our 
results, we assess the mass of AC dust that can grow from the 
C6H6 and C2H2 masses in the radius range of interest, compare 
the derived AC mass to the gas mass, and derive dust-to-gas mass 
ratios. It has been found that IRC+10216 is heavily dust-rich 
and has a dust-to-gas mass ratio comprised between 1.4 x 10""' 
and 4 x lO""' based on various studies (e.g., Ivezic & Elitzur 
1996, Groenewegen 1998, Decin et al. 2010a). In Tabled we 
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Table 5. Derived dust-to-gas mass ratio as a function of radius in the dust formation zone of IRC+10216. We list the pre-shock gas 
number density, acetylene and benzene abundances, the total number of C2H2 and C(,H(, molecules, the CeHs yield (= number of C 
atoms initially in C2H2 locked into CeHe), the mass of coronene C24H12, the total number of C2H2 molecules depleted in the growth 
from C(,H(, to C24H12, the corresponding number fraction of C2H2 depleted in the growth from C(,li(, to C24H12 in %, and the mass 
of AC dust formed (see text for detail). The region of efficient AC dust formation corresponds to numbers in boldface. 
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nCgas)" 


3.65(13) 
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5.76(50) 
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1.55(47) 


YCCfiHs) 








0.2% 


2.7 % 


3.38 % 
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4.89(23) 
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3.24(27) 


3.95(26) 


7.73(25) 
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1.40(48) 
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0% 


0.0003 % 


0.58 % 


6.83 % 


10.16 % 
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0.47 % 


m(AC dust)" 
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4.89(23) 
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1.30 X 10-5 


1.22 X 10 " 


1.75 X 10 


3.28 X 10 5 
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Notes. Gas number densities are in cm ^ and masses in g. The estimated total dust-to-gas mass ratio ranges from 1.2 x 10 ^ to 5.8 x 10 ^ 
(see text for more detail). 
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Fig. 8. Abundances (with respect to H2) of hydrocarbon and aro- 
matic species as a function of radius, starting at r, - 1.2 R*. 
Benzene starts forming at r > 2 reaching a peak abundance 
value at 3 R*. Hydrocarbons such as C4H4, C4H2, and propargyl 
C3H3 are also abundant in the dust formation zone. 

list the abundances and total masses of C2H2 and CeHe formed 
as a function of radius. Spherical symmetry for the inner wind 
is assumed to calculate the shell volume and the total number 
of C2H2 and CeHg molecules. We assume that the growth of 
CgHe to coronene, C24H12, proceeds through the H-abstraction- 
C2H2-addition (HACA) mechanism (Frenklach et al. 1984) with 
the addition of nine C2H2 molecules, and that once coronene 
forms, the growth of grains continues via C24H12 dimerisation, 
coalescence, and coagulation (Cherchneff 2011b), We discard 
other growth mechanisms such as the polimerisation of polyynes 
on aromatic radical sites proposed by Krestinin et al. (2000) be- 
cause the abundances of C4H2 are far too low to foster a signifi- 
cant growth of aromatic structures through this channel. Further 
growth via C2H2 addition on the surface of grains may occur 
but is not considered in our simple estimate. Table |5] indicates 
that there exists a region of high dust-formation efficiency be- 
tween 2.5 R* and 4 R*. This region is characterised by high 
gas densities and low enough temperatures (see Table |2]) to ini- 
tiate aromatic growth, as illustrated in Figure [8] Furthermore, 



the concentration of C2H2 is sufficiently high to secure growth 
up to C24H12 without depleting the C2H2 content too much, as 
illustrated in Table |5] where ~ 10 % of the C2H2 species are 
consumed in the growth to C24H12 at 3 R*. 

To estimate the total amount of AC dust mass formed in the 
inner wind and ejected at 5 R*, we consider a simple formalism 
whereby a parcel of gas moves gradually from 2 R* to 5 R* over 
a certain time span. Assuming a microturbulent velocity of be- 
tween 1 and 5 km s ', a range characteristic of the inner wind 
before gas drag and acceleration by dust (Keady et al. 1988), 
and the stellar pulsation period of Table [1] 43 pulsations (and 
shocks) are necessary for that parcel to reach 5 R* assuming 
a microturbulent velocity of 1 km s"', when the pulsation num- 
ber drops to 9 for a microturbulent velocity of 5 km s '. For 
the modest shock velocities and moderate post-shock conditions 
found in the inner wind of IRC+10216, we assume that the AC 
dust is not destroyed in the hot post-shock gas at each shock 
passage, while PAH species are, but reform in the adiabatic ex- 
pansion phase. Therefore, to estimate the total AC dust mass that 
possibly forms in the inner wind and is ejected at 5 R*, we as- 
sume that the total C24H12 mass is converted into AC grain mass 
and we sum up the derived masses over the number of pulsa- 
tions required to reach 5 R*, interpolating mass values from the 
data given in Table |5] We obtain a total dust-to-gas mass ratio 
that spans the range 1.2 x lO^^* - 5.8 x lO"-'. When compared to 
values derived from observations (1 x lO^^* to 4 x lO""*), these 
numbers are satisfactory and point to a specific region in the 
inner wind where carbon dust grows from PAHs and graphene 
sheets. Our simple derivation is based on a 100 % growth effi- 
ciency of CeHe to C24H12 and a 100 % conversion efficiency of 
C24H12 in AC dust grains. These assumptions clearly maximise 
the AC dust mass value at 5 R*. However, owing to the large 
C2H2 reservoir available in this region, an additional growth pro- 
cess that are not considered in the present derivation includes the 
addition of C2H2 molecules at the surface of graphene sheets, a 
mechanism that would add mass to the final carbon dust bud- 
get of the inner wind. Without pointing to a specific value for 
the dust-to-gas mass ratio, our derived ratios closely agree with 
the dust-to-gas mass ratios derived from observations and are in- 
dicative that AC dust formation at these specific radii efficiently 
proceeds in the inner wind of IRC+10216. 
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Fig. 9. Abundances (with respect to H2) of SiC and SiC2, and 
Mg- or Fe-bearing species involved in the formation of metallic 
clusters as a function of radius. Abundances derived by thermo- 
dynamic equilibrium at r, are moved to radius 1 R* for clar- 
ity. Except for SiCi and metal sulphide cluster precursors, abun- 
dances are quite low in the dust formation zone, and pure metal 
clusters do not form from the gas phase. 

4.6. Other dust precursors: Carbides and sulphides 
4.6.1. SiC2 and SiC 

Silicon carbide (SiC) dust has long been observed in the wind of 
carbon stars through its transition at 1 1 fim (Treffers & Cohen 
1974, Speck et al. 1997) and studies on meteorites have con- 
firmed a AGB origin for some of the pre-solar SiC inclusions 
(Zinner 2007). In the laboratory, the synthesis of SiC nanoparti- 
cles is produced using various experimental methods (e.g., laser- 
induced pyrolysis of gas-phase mixtures of silane, SiH4, and hy- 
drocarbons) and SiC grains are observed to form in the tempera- 
ture range 1600-2000 K. These temperatures are encountered in 
the inner wind, and support the hypothesis that SiC forms from 
the gas phase by chemical kinetic processes similar to those ac- 
tive in SiC synthesis in the laboratory. Despite the nucleation 
processes not being fully understood and characterised, a few 
gas phase species have been identified as intermediates in the 
nucleation of SiC particles and include Si, C2, and cyclic SiC2 
(Fantoni et al. 1991). Silicon dicarbide, SiC2, has been detected 
in the inner wind of IRC +102 16 by mm interferometry with 
abundances that range from < 5 x 10"'^ (Gensheimer et al. 1995) 
to 5 X 10"^ (Lucas et al. 1995). Recent observations with HlFl 
onboard Herschel indicate a SiC2 abundance with respect to H2 
of 2 X 10"^ in the inner wind (Cernicharo et al. 2010a). We thus 
assume that the presence of SiC2 in the dust formation zone re- 
flects the nucleation and condensation of SiC grains at high tem- 
peratures and densities and the role of intermediate plays by this 
species. 

There exist no documented reaction rates for SiC2 forma- 
tion and we rely on the isovalence of silicon with carbon to de- 
rive reasonable rates for specific processes. We also consider the 
identified nucleation routes for SiC clusters and the formation of 
SiC2 and (SiC)2 clusters according to Erhart & Albe (2005). The 
main production process for SiC2 in the gas phase is 

SiC + SiC ^ SiC2 + Si, (10) 

where destruction is commanded by the reverse reaction in Eq. 
[TOl and thermal fragmentation. We also assume that two SiC 



molecules react to form the SiC dimer (SiC)2. The results for 
SiC and SiC2 abundances are shown in Figure |9l The shock at 
r, destroys both the SiC and SiC2 initially present in the photo- 
sphere under TE, but SiC reforms in the post-shock gas at 1 .5 
R* and form SiC2 via the reaction in Eq.[TO]and SiC dimers. At 
5 R*, the SiC2 abundance value agrees well with that derived by 
Cernicharo et al. (2010a), indicating that the molecule may be 
regarded as a by-product of the condensation of SiC clusters at 
small radii. The model shows that both SiC2 and (SiC)2 form in 
large amounts as early as 1.5 R*, and well before the aromatic 
formation zone (2.5 R* - 4 R*) discussed in ^4.51 SiC clusters 
thus represent a high temperature condensate population inde- 
pendent of the synthesis of AC dust. We discuss in more detail 
the consequences of this situation for wind acceleration in |5] 

4.6.2. MgS and FeS 

A strong 30 yum emission band has been observed in carbon- 
rich evolved stars at various stages of their evolution, including 
AGB and post- AGB stars, and planetary nebulae. Observation of 
the band was also reported in lRC+10216 and ascribed to solid 
magnesium sulphide, MgS, for which a low radiative tempera- 
ture comprised between lOOK and 450 K was derived (Goebel & 
Moseley 1985, Szczerba et al. 1999, Hony et al. 2002a). A band 
at 23 fim was observed in two carbon-rich planetary nebulae and 
FeS in the form of troilite was proposed as a possible carrier 
(Hony et al. 2002b). FeS is also responsible for the 23 /im band 
detected in proto-planetary discs (Keller et al. 2002). However, 
there is no observational evidence that this band is present in the 
spectral energy distribution of C-rich evolved stars in general, 
including lRC+10216. Chemical models assuming TE in the in- 
ner wind predict that both MgS and FeS condense in carbon-rich 
environments (Lattimer et al. 1978, Lodders & Fegley 1999). 
Keeping in mind that dust formation is not an equilibrium pro- 
cess in stellar outflows, we test the formation of MgS and FeS 
molecules as the initial gas-phase precursors of MgS and FeS 
grains in the inner wind. We also consider the formation of Mg- 
and Fe- bearing molecular species and clusters. They include the 
pure iron and magnesium small clusters, Fe2 and Mg2 respec- 
tively, although there exists no observational evidence of pure 
metal clusters in AGB environments. However, pure iron grains 
are often proposed as a dust component of O-rich AGB winds 
to account for the required near-lR opacity necessary to acceler- 
ate the wind (Woitke 2006). We also consider the hydrides MgH 
and FeH, and the gas-phase precursors to metal oxides MgO and 
FeO. 

Atomic magnesium is an alkaline earth metal that primar- 
ily reacts with oxygen compounds (e.g., H2O, O2, NO, O3) 
to form magnesium oxides. Although reactions with sulphur- 
bearing compounds are not documented, we expect similar types 
of reaction between Mg and SO to those between Mg and O2, 
owing to the isovalence of sulphur with oxygen. We therefore 
assume that Mg will react with SO according to the reaction 

Mg + SO^MgS + O (11) 

Using again the isovalence of sulphur with oxygen, we assume 
that a reaction similar to the reaction in Eq.[TT] triggers the for- 
mation of molecular FeS. For the formation of iron monoxide, 
FeO, and magnesium monoxide, MgO, reactions of atomic Fe 
and Mg with O2 are considered. These processes have been ex- 
tensively studied (e.g., Akhmadov et al. 1988) and their rates 
well-documented. For both sulphides, we also consider the fol- 
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lowing radiative association reaction as a possible production 
channel 



X + S ^ XS + hv, 



(12) 



where X = Mg or Fe. Studies by Kimura et al. (2005a, 2005b) 
that explore various formation routes to MgS and FeS in the lab- 
oratory show very efficient synthesis from the reaction of the two 
Mg (Fe) and S gaseous phases in gas flash evaporation methods, 
and support the occuiTence of the reaction in Eq.[T2]and its ter- 
molecular analogue process. The formation of Fe2 is described 
by the reaction 



Fe + Fe + M ^ Fe2 + M, 



(13) 



where M is the gas collider A rate for the reaction in Eq. [T3] 
was derived by Giesen et al. (2003) in their study of pure iron 
cluster formation at high temperatures. The radiative association 
reaction between two Fe atoms is also included. Similar types of 
processes and rates are considered for the synthesis of MgT. We 
assume that the reverse reactions of all chemical pathways are 
the only destruction processes operating on Fe2 and Mg2. 

Abundances with respect to H2 for these species are shown in 
Figure |9] except for FeO and MgO which have negligible abun- 
dances in the inner wind (x(MgO, FeO) ~ 1 x 10"'^). Apart from 
MgS and FeS, Mg- and Fe-bearing species have very low abun- 
dances. Specifically, Fe and Mg do not form pure metal clusters 
as Fe2 and Mg2 have low abundances in the dust formation zone. 
According to our model, most of the Mg and Fe initially present 
in the photosphere and at r, stays in atomic form. A moderate 
amount of Mg and Fe is first included in the hydrides MgH and 
FeH (see Figure [3]l but are preferentially included into MgS and 
FeS at r > 2 R* from the reaction of Mg and Fe with SO fol- 
lowing the reaction n Eq. [TT] Similar to H2O, sulphur monoxide 
formation is induced by the release of oxygen atoms in the col- 
lisional dissociation of CO in the hot post-shock gas. 

An upper limit to the total mass of solid MgS and FeS pro- 
duced is derived assuming that all gas-phase MgS and FeS are 
depleted into clusters and grains once the gas has reached the 
temperature regime (T < 450 K) derived by Goebel & Moseley 
(1985). The mass limit for both MgS and FeS is at most ~ 0.2% 
of the carbon dust mass formed between r , and 5 R* . When mod- 
elling the spectral energy distribution of lRC+10216, Ivezic & 
Elitzur (1996) derive a dust composition where MgS accounts 
for less than 10 % by mass. For the carbon-rich post-AGB star 
HD 56126, Hony et al. (2003) found that a MgS mass of ~ 2 
% of the carbon dust mass is necessary to account for the 30 
yum band flux. One would thus expect the MgS mass to represent 
at most a few percent of the AC dust mass in the inner wind of 
carbon stars. This required mass is higher than our upper limit 
by a factor of ten. This discrepancy may be the result of several 
uncertainties in the sulphur chemistry or may point to either a 
nucleation process for MgS that does not occur in the gas phase 
or to a different carrier for the 30 fim band. These various aspects 
are discussed in ^ 

Gas-phase FeS follows an abundance trend similar to MgS, 
as seen in Figure |9] because we assumed that the FeS and MgS 
chemistries were alike. Both Fe and Mg are mainly in atomic 
form and have very similar abundances at TE in the photosphere. 
Both solid MgS and FeS possess identical clustering structures 
going from (XS)2 (X - Mg, Fe) with its planar rhombic struc- 
ture to (XS)4, which quickly reaches a distorted cubic structure. 
However, MgS clusters are unstable in O-rich environments, 
contrarily to FeS clusters. It was proposed by Begemann et al. 
(1994) that composite solid sulphides including both Mg and Fe 
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Fig. 10. H2O abundances (with respect to total gas) as a func- 
tion of radius and pulsation phase (legend). Phases 0.8 and 1 are 
indistinguishable. The grey shaded area corresponds to the zone 
where benzene forms and aromatics grow to form carbon dust 
(see §1431). 



could satisfactorily reproduce the 30 fim band in lRC+10216. 
In particular, a magnesium-iron sulphides whose composition 
ranged from Mgo gFeo.iS to Mgo.5Feo.5S provided the closest 
matches to the band. Such a carrier could indeed be synthesized 
in the dust formation zone of IRC -1-10216 in view of the Fe- 
and Mg-bearing species that form, and the presence of both gas- 
phase MgS and FeS. 

As pointed out before, most of the Mg and Fe in the inner 
wind region is in atomic form. High resolution observations of 
optical absorption lines of several metals by Mauron & Huggins 
(2010) for IRC+10216 also appear to detect atomic metals in the 
gas phase, in direct contradiction with TE condensation mod- 
els for which all metals are depleted in a solid phase once the 
condensation temperature of the solid is reached in the wind. 
However, the derived abundances for iron and calcium atoms 
highlighted some degree of depletion relative to the solar abun- 
dance values. According to this study, the depletion cannot result 
from the trapping in a molecular phase, as corroborated by the 
present results where the abundances of metal-bearing species 
are always less than 1 x lO"*". The partial depletion of iron and 
calcium must thus result from either the incorporation of free- 
flying Fe and Ca atoms in the inner wind during the AC dust 
condensation process at r < 5, or the adsorption of these atoms 
on the surface of dust grains at the lower gas temperatures en- 
countered at larger envelope radii. 

4. 7. Line variability w/fh time 

Our model predicts that specific molecular abundances have a 
strong time-dependence. Some molecular species indeed show 
a strong variability in their abundances as a function of time or 
phase of the pulsation period. The abundance of water is a good 
example of such a variability and its abundance with respect to 
H2 as a function of radius and pulsation phase is shown in Figure 
[TOl At 2 R*, H2O abundances span almost six orders of magni- 
tude over one pulsation period (P = 650 days), and this variation 
should be reflected in the intensity of its high-energy transitions. 
In these deep layers, the transitions are mainly pumped by IR ra- 
diation. Apart from the intrinsic variability of the stellar flux with 
pulsation phase that affects all molecules, the large variations 
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in abundances should have some impact on the high excitation- 
line fluxes. These changes in abundances are the consequence of 
the post-shock non-equilibrium chemistry and the destruction of 
molecular species in the hot post-shock gas at early phases. 

Destruction is more or less severe depending on the species 
and not all molecules behave like H2O. For example, with its 
strong molecular bond, CO does not experience such variations 
and shows a rather constant abundance distribution with radius 
and pulsation phase. Some species (e.g., SiO), show a time vari- 
ation in abundance that increases in amplitude from large radii 
to small radii, but where the variation amplitudes do not span 
more than one order of magnitude. Finally, other molecules (e.g., 
SiCa) show average amplitude variations (of three orders of mag- 
nitude) at small radii but reach time-independent abundances 
deep inside the inner wind (~ 2.4 R*). Although the impact of 
time-varying abundances on the line fluxes is difficult to quan- 
tify without a proper radiative transfer model, we would expect 
some species to show little line variability (e.g. CO) but oth- 
ers to be prone to moderate variability (e.g., SiO) and undergo 
large line-flux changes with time (e.g., H2O). Clearly water is 
an excellent tracer of time variability of high-energy molecular 
transitions and as such, an excellent indicator of shock activity 
and shock-induced chemistry in the dust formation zone. 

5. Summary and discussion 

We have modelled the inner wind of the carbon star lRC+10216 
assuming the periodic propagation of pulsation-driven shocks 
between 1 R* and 5 R* and considering a complete gas-phase 
chemistry that encompasses several chemical families. These 
shocks trigger a non-equilibrium chemistry in the hot post-shock 
gas that leads to the formation of molecules and dust precursors. 
The study points to the following new results and trends appli- 
cable in general to carbon stars: 

- The model confirms the presence of a group of molecules, 
namely CO, HCN, SiO, CS, and H2O, that efficiently form 
between 1 R* and 5 R*. These species are expected to be 
present, albeit in different quantities, in the inner wind of all 
AGB stars, regardless of their C/O ratio, as already proposed 
by Cherchneff (2006, 2011a). The derived abundance values 
agree well with available observations. In particular, the dis- 
sociation of CO by collisions in the immediate post-shock 
gas triggers the formation of atomic O, OH, SiO, and H2O. 

- We have found that some hydrides form a new category of 
abundant and stable species in the inner wind. These include 
AlH, HCl, and HF. Other hydrides do form in large amounts 
at r s but are rapidly converted in the dust formation zone into 
chlorine-bearing species, leaving the wind acceleration re- 
gion with low abundances. The most abundant hydrides will 
be released to the intermediate envelope and participate in an 
active chemistry at larger radii. 

- Once formed from the reaction of H2 with CI, HCl is the 
production agent of other Cl-bearing species such as AlCl, 
NaCl, or KCl. The formation chemistry of chlorine-bearing 
species is thus independent of the C/O ratio of the stel- 
lar photosphere, and Cl-bearing molecules including HCl, 
NaCl, and KCl should also be present in O-rich AGB and 
supergiant stars, albeit with different abundances. AlCl is ex- 
pected to have a lower abundance in the dust formation zone 
of O-rich sources because of the Al depletion in gas- or solid- 
phase metal oxide species. 

- There exists a specific zone that extends from 2.5 R* to 
4 R*, where the closure of CeHg ring occurs through the 



recombination of two propargyl radicals, with a CgHe abun- 
dance peaking at 3 R*. The available C2H2 abundances in 
the gas phase are high enough to secure growth to PAHs 
such as C24H12. These large PAHs will consecutively co- 
alesce and coagulate to form AC dust. The estimated total 
dust-to-gas mass ratio spans the value range 1.2 - 5.8 x 10"^ 
and closely agrees with existing values derived from obser- 
vations of lRC+10216. Within this zone, SiC2 molecules ef- 
ficiently form as by-products of the synthesis of SiC clus- 
ters. Some MgS and FeS species are also produced in the 
gas phase but their abundances are too low to account for the 
26-30 yum emission band. 
- The shock-induced scenario predicts a time-variability for 
some molecular abundances over a pulsation period (e.g. 
H2O and SiO) that should induce a time variability in their 
high excitation line fluxes. Other species show either negligi- 
ble changes in abundances (e.g., CO), or small changes that 
do not affect the molecular line intensity (e.g., SiC2). This 
predicted time variability is the direct result of the destruc- 
tion of species in the hot shocked gas layers. Observations 
of the high-energy transitions of these specific species at dif- 
ferent epochs of the pulsation period would help us to con- 
firm the predicted time variability, and assess the impact of 
shocks on the gas chemistry. 

As reported in § 14.51 we have found that the formation of 
PAH molecules, and both their coalescence and growth to AC 
dust take place in a specific radius range. The growth of ben- 
zene, CeHe, to coronene, C24H12, via the HACA mechanism is 
expected to consume a large part of the benzene rings synthe- 
sised at these radii. At radii larger than 4 R*, inspection of Table 
|4]shows that some benzene rings still form and can grow to larger 
aromatic species as the C2H2 reservoir is still large. However, the 
lower gas densities and temperatures should hinder the coales- 
cence of PAHs to form large graphene structures. A population 
of free-flying PAHs not incorporated in AC dust should thus be 
expelled to larger radii once the wind is fully acccelerated. The 
so-called unidentified infrared bands are observed in carbon stars 
that are part of binary systems (Speck & Barlow 1997, Boersma 
et al. 2006). For the carbon star TU Tau, the UV radiation field 
of the blue companion could excite the aromatics present in the 
carbon star wind. These excited PAHs might include the free aro- 
matics that are synthesised beyond the aromatic formation zone 
highlighted in this study. 

We have also found that SiC dimers form at ~ 1.5 R*, i.e., 
far smaller radii than the aromatic growth region, which implies 
that there is an independent population of SiC clusters at small 
radii. Owing to the extinction properties of SiC dust and a sig- 
nificant decrease in the Planck mean of its extinction efficiency 
for temperatures corresponding to the effective temperatures of 
AGB stars, the dust experiences an inverse greenhouse effect 
for a radiation field characteristic of carbon stars (Oilman 1974, 
McCabe 1982, Yasuda & Kozasa 201 1). Since the pressure force 
acting on dust grains is directly proportional to the Planck mean 
of the extinction efficiency, most of the acceleration of the wind 
is provided by AC dust grains (Cherchneff et al. 1991). On the 
other hand, the inner SiC cluster population should experience a 
minor radiation pressure force and lags behind the AC clusters 
when expelled to larger radii. This situation may be reflected in 
the results of studies of meteorites. Pre-solar SiC grains bearing 
the isotopic fingerprint of the AGB s-process are not included 
in graphite spherules that have an AGB origin but form a sep- 
arate pre-solar grain population (Hynes et al. 2007). This iso- 
lation of the SiC presolar grains may directly result from the 
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non-equilibrium chemistry in the post-shock gas that produces 
two main dust populations, namely SiC and AC grains, at two 
distinct positions in the wind acceleration zone. 

No firm conclusions about the production of MgS or Mg-FeS 
dust grains in IRC+10216 can be drawn from the present results, 
as they can be interpreted in many ways. Firstly, the chemical 
model may underestimate the MgS and FeS abundances by a 
factor of ten or more because too much atomic sulphur is trapped 
in SiS, as indicated by the slightly higher SiS abundances listed 
in Table [3] relative to those derived from observations. A small 
amount of S not locked in SiS would result in higher SO abun- 
dances and in a larger amount of MgS and FeS owing to the 
large reservoirs of free atomic Mg and Fe in the inner wind. 
The MgS and FeS clusters would then be produced from the gas 
phase at r > 2 R*, and the resulting estimated MgS dust mass 
could reach the few percent of AC dust mass necessary to ex- 
plain the 30 fim band. That the present model forms gas-phase 
FeS and MgS with similar efficiencies indicate that it may form 
composite Mg-Fe sulphide clusters instead of pure MgS clusters, 
as proposed by Begemann et al. (1994). According to Kimura et 
al. (2005a, 2005b), a MgS (FeS) nucleation from the gas phase 
produces spherical cubic MgS (FeS) clusters instead of the elon- 
gated network-like grains that form when gas-surface reactions 
are involved in the nucleation process. Hony et al (2002b) stud- 
ied the effect of dust shape and temperature on the band shift 
in wavelength with a shift towards 26 yum when spherical and 
hot grains were considered. The emission band in IRC+10216 
clearly peaks around 27 fim in the ISO spectrum, pointing to 
a possible formation pathway from gas-phase chemistry, as we 
have described in the present study. Secondly, if the MgS (FeS) 
abundances are indeed low in the dust formation zone, they point 
to: 1) a synthesis mechanism for MgS or Mg-FeS involving gas- 
surface processes, or 2) alternative carriers for the 30 fim band. 
However, if MgS grain formation occurs at lower temperatures 
and on the surface of already produced dust grains (see the com- 
prehensive studies by Men'shchikov et al. (2001) and Zukhovska 
& Gail 2008), their growth is hindered by the lack of available 
atomic sulphur, which is chiefly depleted in SiS and CS in the 
dust production zone. Therefore, a formation scenario involving 
surface chemistry would also require a mechanism that could re- 
turn sulphur to the gas phase just after the acceleration of the 
outflow. Another explanation is that MgS is not the carrier of the 
30 jum band. A critical assessment of all previous MgS studies 
was made by Zhang et al. (2009) who pointed out that the mass 
of MgS derived from the emission at 30 fim violated the avail- 
able abundances of Mg and S in the stellar atmospheres, owing 
to the use of improper optical constants for MgS in the optical 
and UV wavelength domains. Alternative solids have been pro- 
posed (e.g., hydrogenated amorphous carbon, HAC, Grishko et 
al. 2001). A fresh reinvestigation of the carrier of the band in 
IRC+10216 would be extremely useful coupled to observations 
of high energy transitions of sulphur-bearing species to constrain 
the sulphur reservoir in the dust formation zone. 

Finally, the results presented in this study are not unique 
to IRC+10216 and similar trends should apply to other carbon 
stars as well. For the specific case of water, H2O has now been 
detected in several carbon stars (Neufeld et al. 2011b) where 
formation processes similar to those described in Cherchneff 
(2011a) and in the present study take place. However, the H2O 
abundance certainly varies from source to source depending on 
the various parameters that are entangled with its formation in 
a complex way. For example, the shock strength affecs both the 
destruction of molecules and the creation of free atomic oxy- 
gen. Therefore, one would expect lower water abundances to be 



created by a mild shock than a strong shock, but the opposite ac- 
tually occurs. We have modelled the chemistry induced by a 10 
km s"' shock at r , and compared our results with those of the 20 
km s ' shock, finding that more SiO and H2O molecules were 
produced. Because less CO is destroyed by a mild shock, the 
SiO formation depends on the reaction of Si with CO. Therefore, 
the OH radical is free to form H2O, and combined with the less 
efficient destruction of molecules in the post-shock gas, more 
water is formed. Hence, the many parameters affecting the for- 
mation of water in carbon stars include the shock strength as 
well as the chemical composition, both the photospheric gas den- 
sity and temperature, the gas-phase chemistry of the Si and S 
chemical families, and the amount and type of dust that forms. 
Water abundances is thus expected to vary greatly in carbon stars 
despite its synthesis in these objects having been proven both 
observationally and theoretically. A similar conclusion may be 
drawn for other species such as SiO. Combined observations of 
several high excitation transitions of H2O and SiO molecules 
would be very instructive in this regard to more clearly under- 
stand the chemical processes responsible for the formation of 
O-bearing species and characterise the water content of carbon 
stars on a global scale. 
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Table .1. The chemical reaction network ordered by type of reactions. The reaction rate is given in the Arrhenius form A x (T/300)" 
X exp(- Efl/T), where T is the temperature in Kelvin. The related parameters are given in each column as follows: Column 1 - reaction 
number; Column 2 - reactants and products; Column 3 - A parameter (in cm-' s"', cm*" s"', and s"' for bimolecular, trimolecular, 
and unimolecular processes, respectively); Column 4 - n parameter; Column 5 - activation barrier E„ in Kelvin: Column 6 - rate 
reference. 



TERMOLECULAR 

3B1 H + H + H, — ^ H2 + H2 8.85x10-" -0.6 NIST 

3B2 H + H + H — ^ H2 + H 8.82 x 10^" NIST 

3B3 H + H + He Hj + He 4.96 x 10^'^ NIST 

3B4 H + O + M OH + M 4.36 x 10"" -1.0 NIST 
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